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Abstract 

Recently DELPHI Collaboration reported new data on Bose-Einstein correla- 
tions (BEC) measured in e + e~ — ► W + W~ events. Apparently no enhancement has 
been observed. We have analyzed these data including final state interactions (FSI) 
of both Coulomb and strong (s-wave) origin and found that there is enhancement 
in BEC but it is overshadowed by the FSI which are extremely important for those 
events. We have found the following values for the size of the interaction range (3 
and the degree of coherence A: /3=0.87 ±0.31 fm and A=1.19 ±0.48, respectively. 

1. Introduction: Recently DELPHI Collaboration has reported an interesting result 
on the Bose-Einstein correlations (BEC) in e + e~ — > W + W~ events at -y/s = 172GeV Q. 
Namely, analyzing their data by the following formula 

n ±±/ n bg(+~) = c [ x + Aexp(-/3 2 Q 2 /2) ], (1) 

where the denominator N BG ^ + ~^ denotes the contribution of the different charge pair, (c 
and A are the normalization factor and the degree of coherence, respectively), authors of 
Ref. found that 

c = 1.0 , A = -0.22 ± 0.2(statics), and (3 = 0.5 fm, 

i.e., apparently there is no enhancement due to the Bose-Einstein effect in e + e~ — > W + W~ 
events, most probably because measured pions come from different W's (i.e., different 
" sources" ) . 

On the other hand, different approach (weight function method) to the BEC in e + e~ — > 
W + W~ events used recently in Ref. lead to slightly different results, namely A was es- 
timated to be equal 0.032 ~ 0.146Q (See Table I). 
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1 Another scheme (shift of the final state momentum) is recently argued in Ref.0 in order to also 
determine the VK-mass. 



Because in these reactions the space-time distances between W's are very small (of 
the order 0.5 fm) the usual Coulombic and strong final state interactions(FSI) between 
7r ± -7T ± pairs should be specially important here and should be taken into account in order 
to estimate correctly and precisely (3 and A. Therefore, we would like to re-analyze data 
of |l| by formula containing the FSI (of both the Coulombic and strong origin, the later 
in s-wave) proposed by us in Ref. ||. 

In the next section, a theoretical formula for correlation of n^-it^ pairs are derived. In 
Sec. 3 we derive also a theoretical formula for 7r + -7r~ pairs. Using these two formulas we 
analyze in Sec. 4 e + e~ — > W + W~ events showing clearly the roles of strong and Coulomb 
final state interactions in the BEC. We also predict in this Section correlation functions 
for N ±ziz /N BG and N + ~ /N BG . Our concluding remarks are presented in the last Section. 



2. A theoretical formula for n ± n ± pairs including the FSI: 

To describe correlations of pair of identical bosons we need to symmetrize the correspond- 
ing total production amplitude, which reads now: 

A 12 = 4|[*c(k, r) + *g(k, r) + $ st (k, r) + $f t (k, r)], (2) 

where superscript S denotes the symmetrization of the wave functions. The functions 
\I / c(k, r) and $ st (k, r) stand for the Coulomb wave function and the wave function in- 
duced by the strong interactions, respectively. 

The Coulomb wave function of the identical ir-ir scattering with momenta p\ and p 2 
is given as || || : 

*c(k,r) = r(l + ^)e- m?/2 e ik - r F(-zr/;l;iA;r(l-cos0)), (3) 

where F denotes the confluent hypergeometric function, 2k = Q = pi — p 2 , and 
r) = m, n a/2k. For calculations of the confluent hypergeometric function F in Eq. (§), 
there are two methods. One is to use formulas given in Ref.0. The other is to use the 
subroutine program 'WWHITM ' in the CERN program library (See Refs.[fj], ||). 

The strong interaction in the Coulomb field is given by 

$ st (k, r) = / ° WeXP[ ^ r r ' r/ln(2fcr))] , (4) 

/o(0) = ^exp[2zar g r(l + n7)][ exp(2td^) - 1 ] . 



Here yG(r]) = y ^^p^p ~[ is the square root of the Gamow factor which is assumed to be 
a normalization factor of the strong wave function^ and 5^ is the the phase shift for the 
2 This normalization factor is introduced to remove divergence of fo(0) at k = 0. 



7T- 



7T (1=2) interaction parametrized as[][' 



(2) 



2 1.0 + 0.5Q 2 



(5) 



We obtain therefore 

A r±± j/]\[ BG ( mix ) 



p(r)d 3 r\A 12 \ 2 , 
G(rj) [ (1 + A 1C ) + (E 2B + A EC ) + J Cst + I st 



(6) 



where 



(1 + A 1C ) + (E 2B + A EC ) = 
E 2B = J d 3 rp(r) exp[-zQ • r 
/st = GCk) I rf3r P( r )|^|[ <l> st(k, r) + $f t (k, r 



/ d 3 rp(r) |-^ c (k, r) + *g(k, r '* 



'Cst 



G(kY 



Re d 3 rp(r) {* c (k, r) + *g(k, r)} x {$ st (k, r) + $f t (k, r)}* ,(7) 



where N BG ( mix ) is number of 7r + 7r pairs obtained by event mixing method. In present cal- 



culations we use Gaussian source function, p(r) 
of which is given by: 



1^3 



27I-/3' 



exp(-^r), the Fourier transform 



E 2B = exp(-(3 2 Q 2 /2). 

To obtain the ratio N ±:iz /N BG( - mix ^ the parameter A should be introduced into it in the 
usual way. Notice that one more parameter, the additional normalization factor c, is also 
introduced by hand. Our final formula with these parameters is then given as: 

N ±± /N BG( mi *)( Q = 2f£) = G ( v) x c(1 + + + j cst + /st) 

E 2B 



X 



1 + A 



1 + Aic + A EC + Ickt + 1st 



It should be noted that the normalization parameter c and an effective degree of co- 
herence, i.e., the denominator of the ratio E 2B /(1 + Aic + Aec + ^Cst + ^st), are related 
to each other. 



3. A theoretical formula for n + ix~ pairs including the FSI: Using the following 
amplitudes of the strong interaction for the n + 7T~ scattering (I =2 and 0) and changing 
the sign of the factor rj = ma/2q, one can estimate correlations for the 7r + 7r~ pairs. The 
strong amplitude for 7r + 7r~ — > 7r + 7r~ scattering is decomposed as 



|/o w (*) + f/o ( >). 



(9) 



3 It is numerically confirmed that calculations for the BEC using Eqs. (|3|) and (fj) is equivalent to 
that using a solution of Schrodinger equation with Coulomb and p-meson exchange potential [H. 



Here, we have assumed that contribution of 7r 7r° — > 7r + 7r~ channels is small. The 
phase shifts 5^ for 1 = which are obtained in Refs. U-|RJ are parameterized by the 
polynomials (cf., Fig. 1). 

5^(Q) = -75.63 Q + 2345. Q 2 - 12588. Q 3 + 23872. Q 4 + 7431.5 Q 5 

-77590. Q 6 + 87705. Q 7 - 30866. Q 8 [deg] (for Q < 0.985 GeV) 
S^iQ) = -53061. + 179950. Q - 227210. Q 2 + 127110. Q 3 

-26558. Q 4 [deg] (for 0.985< Q < 1.350 GeV) 
^ 0) (Q) = +3978.5 - 7988.6 Q + 5608.9 Q 2 

-1269.3 Q 3 [deg] (for 1.350< Q < 1.750 GeV) 
5^(Q) = 388.05 [deg] (for 1.750 GeV< Q ) 

(10) 

Our formula for the 7r + 7r~ pairs is given as 

vl/ total (k,r) = vl> c (k,r) + $( 2 ' 0) (k,r), 

N (+~)/ N BG(m lX ) = J d 3 r p(r) |^ total ( k)r )|2 ; (U) 

It should noted here that Eq. (|TT|) gives a prediction for 7r + 7r _ correlation function caused 
by the FSI. Using now Eqs. (|8]) and (|Tl~D, one obtains the following theoretical formula of 
the BEC with unlike-particle reference [1|] as 

7r ± 7r ± /7r + 7r- = Ar ^ / /ArR _ r . . . 12 



4. Analyses of e + e — > W + W events by means of Eq.( |l2|) : We have now 
analyzed data on e + e~ — > W + W~ events at \fs =172 GeV M using Eq.(|l^). A result 



of minimum-x 2 fit is shown in Fig. 2. Because the scattering length parameter in Eq. (|5]) 
was fixed as a Q = — 1.2OGeV _1 and the normalization parameter c was fixed as c = 1, 
we can estimate the remaining two parameters as: 

• A = 1.19 ±0.48, /3 = 0.87 ± 0.31 fm, ( x7 NDF = 14.3/13 ). 

It should be noted that the value of the A parameter obtained by Eq. (|i~2|) is considerably 
different from one estimated in Ref . |1|, 0] . See Table I. 

In order to expose the cause this difference in a more clear way, we also analyze the 
data by the following two formulas: 



N ±±l N BG{mix) 

7T /tt (strong corr. only) = N+ _ /NBG{mix) 

_/V±± / J\[BG(mix) 

ix /n + (Coulomb corr. only) - 



j\H — j]\[BG{w,ix) 



(13) 

a = 

(14) 

Me) = o 



Results of minimum-x 2 fit to the data by means of Eqs. (13) and (|14]) are summarized 
in Table E. We find that the correction caused by the strong FSI change drastically the 
value of the A parameter. The fitting results by means of Eqs. (13) and (|14]) are also 



shown in Fig. 3. Since the Coulomb correction plays important role only in small Q region 
(Q < 0.2GeV), strong interaction is dominant in the Q region observed by authors of 
Ref.05 (cf. Fig.2 and 3.) 

Finally we also predict the correlation functions, N ±:k /]\r BG ( mix ) anc [ jV + ~ /j\f BG ( mix ) 
in Fig. 4. 



5. Concluding remarks: We have analyzed the data of the BEC in e + e~ — > W + W~ 
events by means of Eq. (|i~2"l) , i.e., including full FSI structure (both Coulombic and strong). 
In this way we have found that there is normal Bose-Einstein effect in the correlation data 
because the degree of coherence parameter A included in Eq. ([12]) comes out to be approx- 
imately unity. However, this expected enhancements is dissolved by the complicated 
structure of the FSI (especially important here due to the small space-time distances in- 
volved in this particular reaction). 
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Figure captions 

Fig.l: The phase shifts of 5^ parameterized by Eq. ([II]). See Refs. [9]- [16]. 

Fig.2: Our result of analysis for the BEC data by DELPHI Collab.|L| using Eq. ([TJ). 

Fig. 3: Illustration of the role of the Coulomb and strong final state interaction in 
7T ± 7r ± /7r + 7r _ . 

Fig. 4: Our predictions for correlation of 7r =t 7r =l= pairs (solid) and b) 7r + 7r~ pairs (dotted). 



Table captions 
Table I: Effect of FSI on BEC in e+e" -> W + W~ events. 

Table II: Influences of the Coulomb and strong final state interactions for the (5 and A 
parameters. 



Refs. 


P [fm] A 


X7NDF 


FSI 


Eq.(l) 

Event weighting [2] 

Eq.(12) 

7r ± 7r ± /7r + 7r _ 


0.50 -0.11 ± 0.2 
0.322 ~ 0.459 0.032 ~ 0.146 
( depending on weight function ) 

0.86 ± 0.31 1.19 ± 0.48 


~ 14/14 
14.3/13 


no-FSI 
no-FSI 
Coulomb+ strong 



Table I 



Eq. FSI 


P [fm] 


A 


X 2 /NDF 


Eq.(13) Coulomb corr. 
Eq.(14) strong corr. 


0.35 ± 0.17 
0.87 ± 0.32 


-0.15 ± 0.13 
1.16 ± 0.48 


11.5/13 
14.4/13 



Table H 




Fig. 1 
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Fig. 2 
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eq.(14): strong corr. X= 1.16 +/- 0.48 
eq.(13): Coulomb corr. X= -0.15 +/- 0.13 



0.25 

0.00 0.50 1.00 1.50 2.00 2.50 3.00 

Q [GeV] 



Fig. 3 
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Fig. 4 



